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Abstract 

A simple approach to describe the levitation force measurements 
on melt-processed HTS bulks was developed. A couple of methods to 
determine the critical current density J c were introduced. The aver- 
aged aft-plane J c values for the field parallel to this plane were deter- 
mined. The first and second levitation force hysteresis loops calculated 
with these J c values coincide remarkably well with the experimental 
data. 



1 Introduction 

Superconducting systems with magnetic levitation have long been known and 
the discovery of high temperature superconductors (HTS) highly stimulated 
their investigation but a real interest in them for large scale applications 
appears only with the successful development of the melt-processed (MP) 
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technology PJ. The use of MP HTS in large scale systems such as flywheels 
for energy storage, electric motors and generators, permanent magnets, etc. 
is the most promising HTS application now 0. In this applied region the lev- 
itation force measurements can be considered in two roles: as an information 
source to know more about levitation systems and as a quick technique to 
test HTS samples 0. In many earlier works j|, |5| it has been shown that the 
forces between a PM and HTS sample are closely related with HTS magne- 
tization curves. Vertical levitation force versus vertical distance F z (z) is the 
nearest analog to M(H) dependencies with their major and minor hysteresis 
loops but the complexity of a field configuration in such large scale PM-HTS 
systems makes it very difficult to directly correlate them in general case. The 
problem can be solved by numerical approaches and some of them have been 
successfully used ||. The numerical approaches are undoubtedly useful to 
evaluate the real system parameters but usually need too much computer 
resources to be applicable to direct HTS sample investigation. To perform 
such an investigation an analytical evaluation is more wished for. 

Two limiting cases of HTS structure have been considered as analyt- 
ical to calculate the dynamic parameters of an idealized system, a point 
magnetic dipole over an infinite flat superconductor. The first one is the 
case of 'granular superconductor' which can be modeled by a set of small 
isolated superconducting grains |7|]. It was shown the usual granular HTS 
obtained by standard sintering techniques can be described very well within 
this model [0, |8j. The second one is the case of an 'ideally hard supercon- 
ductor' H [R], [11], 12 1 . It was shown recently || 10, 11] that dynamics in a 



wide variety of levitation systems can be described in terms of surface screen- 
ing currents which screen alternating magnetic field component due to PM 
displacements like it would be for an ideally hard superconductor, a super- 
conductor with infinite pinning forces. For an infinite flat superconductor 
the frozen-image method was introduced as an illustration of simple ana- 
lytical calculation of forces acting in the system. A perfect agreements with 
experiments were found by us for small PM resonance oscillations frequencies 
|IIJ and, recently, by Hull and Cansiz |TB|] , for both vertical and lateral force 
components. 
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2 Approach description 



The feasibility of the 'ideally hard superconductor' approach is that the pen- 
etration depth S of alternating magnetic field is much less than system di- 
mensions [K| . To calculate the stiffness or resonance frequencies the limit 
5 — > can be used, but it was shown that taking into account the finite values 
of S it is possible to calculate the a.c. loss [|11]] and even recover critical cur- 
rent density profiles within 5 depth from a.c. loss measurements [Q. In this 
paper we present such an approach of levitation force calculation (including 
its hysteresis behavior) for superconductor with finite values of critical cur- 
rent density J c and simple methods to obtain J c values from levitation force 
experimental data. 

A PM placed over ideally hard superconductor induces at its surface 
the screening currents j = (c/47r)n x b r , where n is the surface normal and 
b r is the tangential magnetic field component at the surface (the normal 
component b n at the surface is zero). From the symmetry, for an infinite flat 
surface WO 



h r = 2h ar , (1) 

where b a is the variation of the PM magnetic field B a due to its displacements 
in respect to initial field cooled (FC) position: b a = B a — B aFC [[Hj]. For the z- 
axial symmetric configuration r = (rsind, rcos8, z) where only jg component 
is induced, for the vertical force acting on PM from the screening currents 
one can write 



id 



rb 2 ar (r)dr. (2) 



This is ideal force which can be readily calculated just from known tan- 
gential component of PM field. The Eq.(Q) is obtained from zero-depth 
screening currents approach that we will call a zero approximation of real 
PM-HTS systems. Within this approximation any configuration of such sys- 



tems can be calculated numerically [14| but to describe hysteresis phenomena 



next-order approximations have to be considered. 
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In the second stage (a first approximation) we will examine a model where: 
(i) S is finite but still much less than system dimensions L, (ii) the critical 
state model is applicable to these samples, and (iii) critical current density 
is constant. The applicability of the critical state model to melt-processed 
HTS has been proven in many experiments |3|, 111, |12| and is quite acceptable 
here. The first condition on 5 can be written as 



*(r) « B ar (r) (^}) '-L, (3) 

and, because of 5 oc b r , can be satisfied anyway limiting the minimum dis- 
tance between PM and HTS surface. One can estimate L rs z + d/2, where 
z is, here and below, the distance between PM and HTS surface and d is the 
PM thickness. 

As for condition on J c , it is far from reality by itself since the critical 
current density usually depends on both magnetic field and space coordinate. 
But as it will be shown below we can accept this condition for levitation force 
measurements. This just means that in the next relation for j, the surface 
density of screening currents, 



j ( r , z) = 5 (r, z) J c = —b ar (r, z) , (4) 

J c can be treated as a coefficient between j and S, a coarse-grained flux pen- 
etration depth averaged over L scale. The j(r) function does not depend 
on field history but only on PM position z in the same way as b ar (r), the 
distribution at the HTS surface of the PM field variation, that after cool- 
ing is a function of r and z. Thus, in the considering approximation, the 
function 5(r, z) formally does not depend on field history but means the flux 
penetration depth at the first PM descent only. 

Next, if we use a protocol of PM motion according to which it moves 
between two points: the initial or FC point z max that is included in b ar (r, z) 
function as a condition b ar (r,z max ) = (z max = oo for ZFC case), and the 
lowest point z m i n , the current distributions in the depth z of superconductor 
are the following. After the PM first stop and beginning to go up (the first 
ascent), the depth of the layer where currents flow remains constant and 
equal to its maximum value 5 max = 5(r, z min ) but there are two regions with 
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opposite currents. The opposite flowing current penetrates from the top at 
the depth 8^ that can be obtained from (4) 



<5 T (r, z, z min ) = -(S(r, z min ) - 8(r, z)). 



(5) 



Its maximum value is 8 



\max 



the second descent there are three regions with 



) = 8{r,z min )/2, so during 
J c for < C < $i , -J c for 



8 l < ( < 8^ m ax , and + J c for <5 Tmaa; < C < $max , where 8 l 



8(r, z)/2 also 



does not depend on z min . If one can neglect flux creep for times greater than 
the descent- ascent time, any other ascents are equal to the first one and any 
other descents are equal to the second one. Any other current distributions 
for other protocols, for example to describe minor hysteresis loops, can also 
readily be obtained within the scheme above. 

Applying this scheme to calculate the vertical forces during the first de- 



scent F(z), the first and the next ascents F^(z, 
next descents F^(z,z min ) one can write 



and the second and the 



oo S(r,z) 

2n f 

F{z) = — J c rdr 

C o n 



d(b ar (r,z + C). 



(6) 



F T (z, z min ) = —J c J rdr J dC, - J dC, 





b ar (r,z + C). 



(7) 



Fi(z,z min ) = — J c I rdr 



2tt 

c 



OO Smax 5 m ax/2 S/2 

J rdr J d(- J d(+ J d( 

5/2 



§max /2 



b ar {r,z + (). (8) 



The functions 8(r, z) depend on J c in according to the above equations 
((H)) © an d below) and for J c — > oo all these forces become equal to Fid(z). 

Remaining within the condition (|3]) we can approximate the integrals 
over z from the formulas (§D-(BP by multiplying the depth of the layer where 
current flows by the field bar in its center. It is easy to show that within the 
above approximation the formula (pi), for example, can be rewritten as 
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F{z) = J rb ar (r,z)b ar [r,z + -\dr, (9) 
o ^ ' 

which highly increases the calculation speed. 

3 Experiment 

To check the applicability of the above consideration to real MP HTS we 
used a standard experimental setup on levitation force measurements [||]. 
The SH1C05 disk shape PM was 15 mm in diameter and 8 mm in thickness 
(the effective thickness with ferromagnetic holder that was evaluated from 
real PM field configuration was 12.7 mm) with averaged axial magnetization 
of 47rM = 9236 G (the field measured by Hall probe in its center at the 
distance of 0.8 mm from its bottom surface was 3350 G). The magnetic 
field of the PM was calculated as field of a coil with the same dimensions 
and with lateral surface current density J = cM. All measured samples 
were melt-processed HTS of 30±0.5 mm in diameter and 17.5±0.5 mm in 
thickness. The distance z between PM bottom surface and HTS top surface 
varied from z max = 400 mm (that can be considered as ZFC case) to its 
minimum value z m i n = 0.5 mm. The minimum step of PM motion was 75 
mm. The accuracy of force detecting was 15 mN. Within this accuracy the 
experimental data were reproducible for every sample. Fig.l represents the 
first and second hysteresis loops (the first and second descent and ascent) for 
two samples. 

Within the above approximation we have only one parameter, J c , for 
forces (pp-(§|) (or (0) and analogous ones) to be fitted to the experimental 
ones F exp (z). To do this, we have to choose one of these functions and one 
point Zi, and solve the equation 



F(J c ,z i ) = F exp (z i ). (10) 

The forces calculated from formulas (§])-(§) with the J c values obtained 
from (|ToD in z min point are also represented in Fig.l by solid lines. The 
forces calculated from the formula (pf) and from analogous ones practically 
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coincide with the above in the F(z) plot scale. A good agreement between 
the experimental and calculated F(z) dependencies demonstrates the above 
approximation is correct. 

Nevertheless, the discrepancy between the experimental and calculated 
forces still exists and is lager than the experimental accuracy. One of the 
most likely reasons is a variation of J c with depth and field. Fig. 2 shows the 
values of J c versus maximum value of B r (z) at the HTS surface for two HTS 
samples. The data were obtained by solving Eq. (|T0|) . Open symbols represent 
the solution for the function (||), and solid symbols represent the solution for 
the function (^). The solid line in Fig. 2 with respect to right axis represents 
the dependence of B rmax {z). For a perfectly uniform sample with c-axis 
exactly perpendicular to the surface such a dependence of J c (B rmax ) would 
be uniquely determined by the dependence of J^ b (B ab ), the critical current 
density flowing in afr-plane versus the magnetic field parallel to this plane. 
But for real melt processed samples, it is more reasonable to assume that 
the dependencies of J c (B rmax ) in Fig. 2 are mostly caused by space variations 
of critical current density. The steep slope of the curves in Fig. 2 at low 
field, which caused the maximum in the apper curve, is reflated with the 
finite diameter of the HTS samples and shows the lower field limit for given 
configuration. 

4 Another simple method 

There is a possibility here to introduce a visual simple method to evaluate 
J c . It is understandable that in spirit of the above consideration a shift Az 
(see Fig.l) of the first descent experimental curve with respect to ideal one 
has to be proportional to an average penetration depth. From the condition 
F id (z + Az) = F(J C , z) and Eq.(|]) one can readily obtain 

c b 

5 « 4Az, or J c « (11) 
8nAz v ' 

The values of J c evaluated in such way are also represented in Fig. 2 
and show a good agreement with ones determined before. The experimen- 
tal error a,j c that is shown here was estimated from the formula <j,j c /J c = 
a f{AF j dz)~ x I Az which assumes the maximum error is caused by the force 
measuring: op « 30 mN. 
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5 Conclusions 



In summary, we have considered the approach, which we call the "first ap- 
proximation", to describe levitation force data. The term "first" implies 
that we consider a case in which such parameters as flux penetration depth 

5 or normal component of magnetic field at HTS surface b n are already not 
zero, as it is for ideally hard superconductor 0, but small enough: S <C L, 
b n b r . Within this condition the methods to calculate J c , the critical cur- 
rent density, which we have introduced in the paper are exact. Remarkably, 
the approach works well even beyond this condition, when 5 ~ L, h n ~ b r . In 
this region the methods become empirical. The J c value that can be obtained 
by the methods is averaged over L scale critical current density in afe-plane 
for field parallel to this plane: J c = (J^ b (B || ab)). L scale depends on the 
size of a magnet we use. 
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7 Figures Captions 



Fig. I. Experimental (symbols) and calculated (solid lines) data on the first 
and second hysteresis loops of the vertical levitation force vs. distance z 
between PM and HTS surface. Dashed line represents the force for an ideal 
superconductor. 

Fig. 2. The values of averaged critical current density versus maximum value 
of B r , the magnetic field tangential component at the HTS surface, obtained 
by different methods. The solid line with respect to right axis represents the 
dependence of B rmax (z). 
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This figure "FIG1.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/9905032vl 



This figure "FIG2.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/9905032vl 



